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Abstract 

Within the approach of Supersymmetric Quantum Mechanics associated with the 
variational method a recipe to construct the superpotential of three dimensional con- 
fined potentials in general is proposed. To illustrate the construction, the energies of 
the Harmonic Oscillator and the Hulthen potential, both confined in three dimensions 
are evaluated. Comparison with the corresponding results of other approximative and 
exact numerical results is presented. 
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1. Introduction 

Super symmetric Quantum Mechanics, (SQM), was conceived in 1981 as the sim- 
plest field theory that would provide understanding of the supersymmetry breaking 
mechanism of higher dimensional field theories, [1]. Since then this formalism has be- 
come a new field of research and a great deal of work, both analytical and numerical, 
has been done to get better knowledge of the exactly solvable, the partially solvable, 
the isospectral, the periodic and the non-exactly solvable potential problems, (see 
[2] for the latest review). In particular, few of these latter have been investigated 
through a novel methodology based on the association of the variational method 
with SQM formalism. The conjunction of this well known quantum mechanism with 
SQM has given a simple framework to investigate spectral problems of Hamiltonian 
systems associated to the three dimensional Hulthen, Morse and Screened Coulomb 
potentials, [3]- [5]. The main basis of the method consists in setting an Ansatz for the 
superpotential. Based on the superalgebra intrinsic of SQM we are able to evaluate 
the wavefunction which naturally contains free parameters introduced by the Ansatz. 
Having in hands our ideal wavefunction, the calculations to be done are the usual vari- 
ational calculations in which the parameters are varied until the energy expectation 
value reaches its minimum. 

More recently the application of this methodology to exploit the energy spectra 
of a confined quantum mechanical system gave encouraging results: the confined 3- 
dimensional hydrogen atom showed results compatible with those obtained from other 
approximative methods and numerical calculation, [6]. 

Motivated by these results and considering the widespread interest in quantum 
confined systems, particularly in the studies of semiconductor heterostructures (nanos- 
tructures like quantum dots, quantum wires, quantum wells), [7]- [9], as well as in field 
theory, [10], in this work we propose a systematic construction of the superpotential 
for 3d confined systems, that yields, through SQM, to the trial wavefunction pursued. 

To illustrate the construction, the technique is applied here to get the Is and 2p 
energy states for the confined Harmonic Oscillator and for the Is, 2p and 3d energy 
states of the confined Hulthen, both in three dimensions. They are very good when 
compared to recent results obtained from other approximative methods as well as 
exact numerical results, [11]- [15]. 

In what follows, in order to be self consistent and to set up the notation, we 
present a brief description of SQM and its association with the variational method. 
At this point we introduce our recipe for the superpotential in the general case. Then 
we show the results of the application of this formalism to the confined Harmonic Os- 
cillator and Hulthen potential. Comparison with the corresponding results of other 
approximative and exact numerical results is presented in the tables, followed by a 
discussion and our conclusions. 

2. Supersymmetric Quantum Mechanics and the Variational Method 

The starting point is the factorization of a Hamiltonian H associated to a given 
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potential V(r), in terms of bosonic operators A ± and the lowest energy eigenvalue 
E . In h — c — 1 units, this is written as 



H = -~ + V(r)=A^A- + E . (1) 

We recall that for the cases we are considering the potential is symmetric so that 
equation (1) is the radial equation and V(r) includes the barrier potential term. The 
bosonic operators are defined in terms of the so called superpotential W(r), 

Thus, substituting the bosonic operators into equation (1) we conclude that, as a 
consequence of the factorization of the Hamiltonian H, the Riccati equation must be 
satisfied, 

W 2 -W' = 2(V(r)-E ), (3) 

where W — Also through the superalgebra, the eigenfunction for the lowest 
state is related to the superpotential W by 

rr 

q? ( r ) = Nexp(- / W(f)df). (4) 
J o 

The above scheme works perfectly well if the potential is exactly solvable. How- 
ever, if the potential is non-exactly solvable, the Hamiltonian is not exactly factoriz- 
able, i.e., there is no superpotential that satisfies equation (3). On the other hand, the 
Hamiltonian can be factorized in terms of an approximated superpotential giving rise 
to an effective potential. This is achieved by making an Ansatz for the superpotential, 
which naturally introduces free parameters in the problem, denoted by the set {//}. 
Thus, through the superalgebra, the wavefunction, (equation (4)), is evaluated. It is 
the trial wavefunction of the usual variational method, denoted by ^ and depending 
on {fi}. Therefore, having in hands the main ingredient of the variational method, 
namely the trial wavefunction, the energy expectation value given by 

can be evaluated and the parameters can be varied until it reaches its minimum value, 
E(p), which is an upper limit of the sought energy level. Indeed, the Ansatz in the 
superpotential drives us to an effective potential V e ff that has a similar form as the 
original potential, i.e., 

V eff (r)= l -(W 2 -W')+E^) (6) 

where W = W(p) is the superpotential for {// = p,}, the set of parameters that 
minimises (5). The above methodology was applied in the search of approximate 
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energy levels of 3-dimensional systems, like the Hulthen, Morse and Screened Coulomb 
potentials, [3]-[5]. 

Concerning the confinement, it is remarkable that for the confined cases we have 
treated so far, the confined Coulomb, [6], and the present cases of the confined Oscil- 
lator and Hulthen potentials, the Ansatz set in the superpotential is restricted to an 
extra confining term added to the superpotential used in the non-confined case, i.e., 

constant 

^confined ''non— confined ^ ~ \ ' ) 

where R is the radius of confinement. The advantage of this choice is that it allows 
us a previous comparative analysis between the original potential (non-confined) and 
the effective potential, containing the infinite (confining) barrier. This is corroborated 
by the fact that in the R — > oo limit the non-confined superpotential is recovered. 
Moreover, the trial wavefunction will naturally vanish at the barrier R = r and not 
as an input, as it is usual in the variational calculation. 

3. The Confined Harmonic Oscillator 

The radial Hamiltonian equation for the Harmonic Oscillator potential, written 
in atomic units is given by 

As the Harmonic Oscillator potential is symmetric, the confinement is introduced 
by an infinite potential barrier at radius r = R. Thus we make the following Ansatz 
for the superpotential 

H/(r) = -^ + -^ + /i 3 r (9) 
r R — r 

which depends of R, the radius of confinement, and of three variational parameters, 
Hi, H2 and h^. As stated in expression (7), the first and the last terms are already 
known from the non-confined case, [2]. The second term deals with the confinement. 
This is explicitly shown through the effective potential. From equation (6), it has the 
following form 

v eff(r) = Y r + /^3-y + 

Mfa - 1) fllp.2 /^3f nn x 

2{R-rf - V(R—7) + lR—^) ^ ( } 

which is clearly infinite at r = R, as expected for a confining system. Notice that 
the effective potential is evaluated for the values of the set of parameters {//} that 
minimise the energy. Our trial wavefunction for the variational method is obtained 
from the superalgebra through equation (4), using the superpotential given by the 
Ansatz made in equation (9). It is given by 

tt„(r) = fi 2 , // 3 , r) ocr w (R - r)^ e^/ 2 . (11) 
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It depends of three free parameters, //i, jj, 2 and jj, 3 and vanishes at r = R. We use 
this trial wavefunction to calculate the energy expection value, given by the equation 
(5), evaluated with the Hamiltonian (8). Its minimisation with respect to the three 
parameters gives E(p,), which, from now on, we call E vsqm . The results are given in 
the tables below for different values of the confining radius R and states Is and 2p, 
corresponding to the values of / = 1 and 2 of the angular momentum, respectively. 
Comparison is made with exact numerical, E exact , perturbative, E pertl [11], other vari- 
ational, E var , [12], WKB, E wkb) [13] and modified WKB results, E centri ,[15], through 
the percentage errors, 

x \E 'exact E VS q m | (-tn\ 

Ovsqm = p A) 

t-'exact 

and 

\E 'exact Epert \ (13) 



Jpert 



E P 



J exact 

The same percetage errors are applied to the other results. We must note that the 
perturbative results lost their accuracy for large values of the radius of confinement, 
for which they are outside the convergence region. For the results of Table 1, R < 1.72 
and of the Table 2, R < 1.84. 

It should be stressed that in the limit of no-confinement, i. e. R — > oo, the vari- 
ational SQM results also agree with the exact non-confined problem, corresponding 
to the removal of the infinite barrier. In this case the energy is exact and is given by 

E = 2n + L + 3/2 . (14) 

Thus for n = and I = 0, state Is, the exact result is -Er^oo — 1-5000, (Table 1), 
and for n = and I = 1, state 2p, the exact result is -Er^oo = 2.5000, (Table 2). 

Table 1. Energy eigenvalues (in Rydbergs) and percentage errors for different values 
of R for the Is state, (/ = 0). 



R 


Eexact 


E 

J -'vsqm 


^vsqm 


Epert 


3pert 


E w kb 


O-wkb 


E va r 


<W 


1.0 


5.0755 


5.0865 


0.22 






5.0627 


0.25 


5.1313 


1.10 


1.5 


2.5050 


2.5104 


0.22 


2.5046 


0.02 


2.5082 


0.13 


2.5265 


0.86 


2.0 


1.7648 


1.7664 


0.09 


1.7588 


0.34 


1.9882 


12.7 


1.7739 


0.52 


2.5 


1.5514 


1.5529 


0.10 






1.5564 


0.32 


1.5567 


0.34 


3.0 


1.5061 


1.5069 


0.05 


1.1532 


23.4 


1.5061 


0.00 


1.5105 


0.29 


4.0 


1.5000 


1.5002 


0.01 


- 4.341728 




1.5000 


0.00 


1.5033 


0.22 


5.0 


1.5000 


1.5000 


0.00 


-48.076319 




1.5000 


0.00 


1.5025 


0.17 


10.0 


1.5000 


1.5000 


0.00 






1.5000 


0.00 
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Table 2. Energy eigenvalues (in Rydbergs) and percentage errors for different values 
of R for the 2p state, (/ = 1). 



R 


Eexact 


E 


dvsqm 




dpert 


^wkb 




E 


Ovar 


E C entri 


Vcentri 


0.3 


112.188 


112.231 


0.04 


112.188 


0.00 


- 




- 


- 




- 


1.0 


10.2822 


10.2847 


0.02 






10.2643 


0.17 


10.3188 


0.36 


10.2876 


0.05 


1.5 


4.9036 


4.9046 


0.02 


4.9034 


0.00 


4.9084 


0.01 


4.9169 


0.27 


4.9068 


0.07 


2.0 


3.2469 


3.2471 


0.01 


3.2434 


0.11 


3.2490 


0.07 


3.2514 


0.14 


3.3081 


1.88 


2.5 


2.6881 


2.6891 


0.04 






2.7079 


0.74 


2.6901 


0.07 


2.6835 


0.17 


3.0 


2.5313 


2.5322 


0.04 


2.3104 


0.73 


2.5310 


0.01 


2.5337 


0.10 


2.5313 


0.00 


4.0 


2.5001 


2.5003 


0.01 


-1.5656 


163 


2.5001 


0.00 


2.5015 


0.06 


2.5001 


0.00 


5.0 


2.5000 


2.5000 


0.00 


-34.0358 




2.5000 


0.00 


2.5012 


0.05 


2.5000 


0.00 


10.0 


2.5000 


2.5000 


0.00 






2.5000 


0.00 











4. The Confined Hulthen Potential 



The radial Hamiltonian equation for the Hulthen potential, written in atomic 
units is given by 

H _ ld2 | *(* + !) fl5 ) 

2dr 2 2r 2 1 - e' Sr K ' 

Once again, as the potential is symmetric, the confinement is introduced by the 
infinite potential barrier at radius r — R. The associated Ansatz for the superpoten- 
tial follows the equation (7) and the expression for the superpotential given for the 
unconfined case, as in [3], 

lv M = ft(lS^) +c ' + ^7 (16 > 

where 

B 1 = -/i 1 (l + l), C 1 = ~ii 1 + J ^. (17) 

This superpotential gives rise through the superalgebra, (4), to the following trial 
wavefunction 

q?^ = i_i 2 ,r) oc (1 - e -" ir )"« e~ Cir {R-rf 2 . (18) 

It depends of two free parameters, fi± and fi 2 and vanishes at r = R, as expected, 
since the effective potential, evaluated by using equations (6), (16) and (17), is infinite 
at this point. Its general form is 

T/ , v fiie-^ 1(1 + 1) pje' 2 ^ a 2 (fl2-l) {l + l)m2e-^ r 



1-e-Pi*- 2 (l-e-^ ir ) 2 2(R-r) 2 (R - r)(l - e~^ r ) 
fiip>2 P>2 1, pi , 1 x 2 



+ n Z n \ + + 7-^) + E(p). (19) 



2(R-r) (l + l)(R-r) 2 y 2 l + l 
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As in the previous case, {p,} is the set of parameters that minimise the energy, given 
by E{p), which we also call, from now on, E vsqm . The results are in the tables 
below, for different values of the parameter 5 and the confining radius R, for the 
states Is, 2p and 3d, corresponding to the values of I — 1,2 and 3 of the angular 
momentum, respectively. Comparison is made with exact numerical, E exact , from the 
1/N approximation, Eiw, [14] and modified WKB results, E centri , [15]. 

For R — > oo exact values of the energy are only for the / = case, and are given 

by 



E 



11 n5 
~2^n ~ ~2 

Thus for 5 = 0.1 and n — 1 the exact result is E R 



(20) 



= -0.45125, (Table 3). 



Table 3. Energy eigenvalues (in Rydbergs) for different values of R, n — 1 and 
1 — 0,1, 2, (states Is, 2p and 3d) and 5 = 0.1. 



R 


1 


Eexact 


-^vsqm 


&vsqm 


E±/N 




E C entri 


&centri 


6.0 





-0.45053 


-0.44945 


0.24 


-0.45109 


0.12 








1 


-0.00865 


-0.00808 


6.59 


-0.00294 


66.0 


-0.00782 


9.60 


7.0 





-0.45111 


-0.45043 


0.15 


-0.45181 


0.16 








1 


-0.04069 


-0.04037 


0.79 


-0.03324 


18.3 


-0.03976 


2.29 


8.0 





-0.45122 


-0.45076 


0.10 


-0.45193 


0.16 








1 


-0.05783 


-0.05762 


0.36 


-0.05293 


8.47 


-0.05510 


4.72 


9.0 





-0.45125 


-0.45090 


0.08 


-0.45188 


0.14 








1 


-0.06728 


-0.06712 


0.24 


-0.06389 


5.04 


-0.06612 


1.72 


10.0 





-0.45125 


-0.45098 


0.06 


-0.45179 


0.12 








1 


-0.07257 


-0.07243 


0.19 


-0.07008 


3.43 


-0.07196 


0.84 


25.0 





-0.45125 


-0.45125 


0.00 


-0.45131 


0.01 








1 


-0.07918 


-0.07915 


0.04 


-0.07920 


0.03 


-0.07921 


0.04 




2 


-0.01390 


-0.01380 


0.72 


-0.01332 


4.17 


-0.01381 


0.65 


50.0 





-0.45125 


-0.45125 


0.00 


-0.45126 


0.00 








1 


-0.07918 


-0.07918 


0.00 


-0.07920 


0.03 


-0.07920 


0.03 




2 


-0.01448 


-0.01448 


0.00 


-0.01450 


0.14 


-0.01450 


0.14 



Table 4. Energy eigenvalues (in Rydbergs) and percentage errors for different values 
of R for the 2p state, (I = 1) and 5 = 0.2. 



R 


exact 


E 

- LJ vsqm 


$vsqm 


E\/N 




Ecentri 


&centri 


8 


-0.01731 


-0.01708 


1.33 


-0.01242 


28.3 


-0.01607 


7.16 


9 


-0.02749 


-0.02731 


0.65 


-0.02428 


11.7 


-0.02612 


4.98 


10 


-0.03339 


-0.03323 


0.48 


-0.03118 


6.62 


-0.03389 


1.50 


25 


-0.04188 


-0.04178 


0.24 


-0.04199 


0.26 


-0.04192 


0.10 


50 


-0.04189 


-0.04189 


0.00 


-0.04196 


0.17 


-0.04191 


0.05 
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5. Discussion and conclusions 

The study of the two 3-dimensional confined quantum systems considered here was 
made through the mechanism based on the association of the variational method with 
SQM, using the recipe given to construct the superpotential, equation (7). For the 
two cases considered, namely, the Harmonic Oscillator and the Hulthen potential, the 
Ansatz, as given by equations (9) and (16), respectively, naturally introduced free pa- 
rameters in the superpotential. Thus, using the superalgebra the wave functions were 
evaluated, showing, of course, explicit dependence on these parameters, equations 
(11) and (18), respectively. They were our optimal wave functions, used to calculate 
the energy expectation values, equation (5), as in the usual variational method. The 
parameters were varied until this expectation value reached its minimum, giving rise 
to the searched energy states. 

The peculiar feature of the approach is that, when doing the Ansatz to the super- 
potential following the recipe of equation (7), we were able to evaluate the effective 
potential, equation (6), and also equations (10) and (19), which resemble the real, 
confining potential under consideration. In both cases, the effective potential is surely 
infinite at the border, r = R, and as a result the wave function, evaluated through 
the superalgebra, equation (4), naturally vanishes at the border because it finds a 
potential barrier that increases until becoming impenetrable at r = R. Furthermore, 
for increasing values of R the border effects vanish and our effective potential be- 
comes the original non-confined potential. The results for the energies improved, 
when compared with the corresponding results of other approximative and exact nu- 
merical results, converging to the results of the non-confined system. On the other 
hand, for small values of R, stronger are the border effects, but even so very small 
deviation from the exact numerical results were found. 

We conclude that our main result is the recipe given for the superpotential when 
confinement takes place, equation (7). The inclusion of the confining term in 
the superpotential was crucial to describe border effects. This prescription showed to 
be appropriate to get good results for the cases considered so far and can be a first 
trial to exploit other confined systems through this method. It is remarkable that the 
approach was applied to different systems. Contrary to the Coulomb potential, [6], 
and the Harmonic Oscillator, the Hulthen potential is not exactly solvable in three 
dimensions. This makes our recipe good enough to work in the general case and 
reinforces the fact that the association of the superalgebra of SQM with the variational 
method gives a simple and useful framework to investigate confined systems in general. 



8 



References 

[1] E. Witten, Nucl. Phys., B 188 (1981) 513 

[2] F. Cooper, A. Khare and U. P. Sukhatme, Supersymmetry in Quantum Mechan- 
ics, World Scientific, Singapure, 2001 

[3] E. Drigo Filho and R. M. Ricotta, Mod. Phys. Lett. A10 (1995) 1613 

[4] E. Drigo Filho and R. M. Ricotta, Phys. Lett. A269 (2000) 269 

[5] E. Drigo Filho and R. M. Ricotta, Mod. Phys. Lett. A15 (2000) 1253 

[6] E. Drigo Filho and R. M. Ricotta, Phys. Lett. A299 (2002) 137 

[7] M. K. Rostov, M. W. Cole and G. D. Mahan, cond-mat/0204094 

[8] F. Gleisberg and W. Wonneberger, cond-mat/0208376 

[9] T. Yildirim and A. B. Harris, cond-mat/0301428 

[10] A. S. de Castro and W. G. Pereira, Phys.Lett. A308 (2003) 131 

[11] V. C. Aguilera-Navarro, J. F. Gomes, A. H. Zimerman and E. Ley Koo, J. Phys. 
A: Math. Gen.16 (1983) 2943 

[12] J. L. Marin and S. A. Cruz, J. Phys. B: At. Mol. Opt. Phys. 24 (1991) 2899 and 
Am. J. Phys. 59 (1991) 931 

[13] A. Sinha, R. Roychoudhury and Y. P. Varshni, Can. J. Phys. 79 (2001) 939 

[14] A. Sinha, R. Roychoudhury and Y. P. Varshni, Can. J. Phys. 78 (2000) 141 

[15] A. Sinha, quant-ph/0205122 



9 



